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ABSTRACT 

We explore the hypothesis that some high- velocity runaway stars attain their peculiar 
velocities in the course of exchange encounters between hard massive binaries and 
a very massive star (either an ordinary 50 — 100 Mq star or a more massive one, 
formed through runaway mergers of ordinary stars in the core of a young massive star 
cluster). In this process, one of the binary components becomes gravitationally bound 
to the very massive star, while the second one is ejected, sometimes with a high speed. 
We performed three-body scattering experiments and found that early B-type stars 
(the progenitors of the majority of neutron stars) can be ejected with velocities of 

> 200 — 400 kms~^ (typical of pulsars), while 3 — 4 Mq stars can attain velocities of 

> 300 — 400 kms^^ (typical of the bound population of halo late B-type stars). We 
also found that the ejected stars can occasionally attain velocities exceeding the Milky 
Ways's escape velocity. 

Key words: Stellar dynamics - methods: N-body simulations - binaries: general 
- stars: neutron - stars: neutron: IRXS J141256. 0+792204 - stars: individual: HD 
271791 



1 INTRODUCTION 

The origin of high-velocity runaway stars can be attributed 
to two basic processes: (i) disruption of a tight massive 
binary following the (asymmetric) supernova explosion of 
one of the binary components (Blaauw 1961; Stone 1991; 
Leonard & Dewey 1993; Iben & Tutukov 1996) and (ii) dy- 
namical three- or four-body encounters in dense stellar sys- 
tems (Poveda, Ruiz & Allen 1967; Aarseth 1974; Gies & 
Bolton 1986; Leonard & Duncan 1990). In the first process, 
the maximum velocity attained by runaway stars depends on 
the magnitude of the kick imparted to the stellar supernova 
remnant [either a neutron star (NS) or a black hole]. For 
reasonable values of this magnitude, the runaway velocity 
does not exceed ~ 200 kms^^ (e.g. Leonard & Dewey 1993; 
Portegies Zwart 2000; see also Gvaramadze 2009). In the 
second process, the ejection velocity could be higher. One of 
the most important and best studied channels for producing 
high- velocity stars is through close encounters between hard 
(Aarseth & Hills 1972; Hills 1975; Heggie 1975) binary stars 
(Mikkola 1983; Leonard & Duncan 1990). Numerical simula- 
tions by Leonard (1991) showed that the maximum velocity 
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that a runaway star (usually the lightest member of the bina- 
ries involved in the interaction) can attain in binary-binary 
encounters is equal to the escape velocity from the surface 
of the most massive star in the binaries. For binaries con- 
taining upper main-sequence stars, the maximum velocity 
of runaways could be as large as ~ 1400 kms~^. The result 
by Leonard (1991) is often invoked to explain the high pecu- 
liar velocities measured (or inferred) for some runaway stars 
(e.g. Heber, Moehler & Groote 1995; Maitzen et al. 1998; 
Tenjes et al. 2001; Ramspeck, Heber, Moehler 2001; Martin 
2006; Gvaramadze 2007; Gvaramadze, Gualandris & Porte- 
gies Zwart 2008, hereafter Paper I; Gvaramadze & Romans 
2008a; Gvaramadze 2009). 

The recent discovery of the so-called hypervelocity stars 
(HVSs; Brown et al. 2005; Edelmann et al 2005; Hirsch et 
al. 2005) - the ordinary stars moving with velocities ex- 
ceeding the Milky Way's escape velocity - attracted at- 
tention to dynamical processes involving the supermassive 
(~ 4 X 10^ Mq) black hole in the Galactic Centre (Gualan- 
dris, Portegies Zwart & Sipior 2005; Baumgardt, Gualan- 
dris & Portegies Zwart 2006; Levin 2006; Sesana, Haardt 
& Madau 2006; Ginsburg & Loeb 2006; Bromley et al. 
2006; Lu, Yu & Lin 2007; Lockmann & Baumgardt 2008; 
Perets 2009). These processes [originally proposed by Hills 
(1988) and Yu & Tremaine (2003)] can result in ejection 
velocities of several 1000 km s~^. Similar processes but act- 
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ing in the cores of young massive star clusters (YMSCs) 
and involving dynamical encounters with intermediate-mass 
(- 100 - 1000 Mq) black holes (IMBHs) were considered in 
Paper I (see also Gvaramadze 2006) to explain the origin of 
extremely high-velocity (> 1000 kms~^) NSs (e.g. Chatter- 
jee et al. 2005; Hui & Becker 2006] and HVSs. Gualandris & 
Portegies Zwart (2007) proposed that exchange encounters 
between hard binaries and an IMBH formed in the core of 
a YMSC in the Large Magellanic Cloud could be responsi- 
ble for the origin of the HVS HE 0437-5439 (cf. Edelmann 
et al 2005; Przybilla et al. 2008a; Bonanos et al. 2008). A 
strong support to the possibility that at least some HVSs 
originate in star clusters rather than in the Galactic Centre 
comes from the proper motion measurements for the HVS 
HD 271791, which constrain the birth place of this early B- 
type star to the outer parts of the Galactic disk (Heber et 
al. 2008; see also Sect. 50 Another example of an extremely 
high-velocity (> 400 kms~^) B-type star originated in the 
Galactic disk is the 5 Mq star HIP 60350 (Maitzen et al. 
1998), whose birthplace lies at about 7 kpc from the Galac- 
tic Centre (Tenjes et al. 2001). 

While the existence of a supermassive black hole in 
the Galactic Centre is widely accepted, no conclusive ev- 
idence has been found for IMBHs. These objects could be 
the descendants of very massive (> 1000 M©) stars (VMSs), 
formed in the cores of YMSCs through a runaway sequence 
of collisions and mergers of ordinary massive stars (e.g. 
Portegies Zwart & McMillan 2002; Portegies Zwart et al. 
2004; see also next Section). Recent numerical studies of the 
evolution of VMSs, however, indicate that these stars can 
lose most of their mass due to copious stellar winds and 
leave behind IMBHs with masses of < 70 Mq (Belkus, Van 
Bever & Vanbeveren 2007; Yungelson et al. 2008), which are 
not large enough to contribute significantly to the produc- 
tion of high- velocity runaway stars (see Paper I) . 

In this paper, we explore the hypothesis (Gvaramadze 
2007) that some high-velocity runaway stars could attain 
their peculiar velocities in the course of strong dynamical 
encounters between hard massive binarie^ and a VMS. In 
this process, one of the binary components becomes gravita- 
tionally bound to the VMS, while the second one is ejected, 
sometimes with a high velocity. Our goal is to check whether 
or not this exchange process can produce early B-type stars 
(the progenitors of the majority of NSs) with velocities of 
^ 200 — 400 kms~^ [ measured for some massive halo stars 
(e.g. Carozzi 1974; Tobin & Kaufmann 1984; Keenan et al. 
1987; Kilkenny & Stone 1988; Ramspeck et al. 2001; Martin 
2006) and typical of pulsars (e.g. Hobbs et al. 2005)] and 
3 — 4 Mq stars with velocities J5 300 — 400 kms~^ (mea- 



^ Note that some authors (e.g. Brown, Geller & Kenyon 2009; 
Tllllch et al. 2009) suggest to use the term "hypervelocity" to 
designate the high- velocity stars ejected solely from the Galactic 
Centre (i.e. via the dynamical processes involving the supermas- 
sive black hole). At present, however, proper motion measure- 
ments are available for only one of the known HVSs (see above) 
so that it is impossible to unambiguously associate the birthplace 
of these objects with the Galactic Centre. In the following we 
will call "hypervelocity" stars all stars with peculiar velocities 
^ 700 kms-l. 

^ We consider a binary as massive if at least one of its components 
is more massive than 8 M.q. 



sured for some late B-type halo stars; e.g. Maitzen et al. 
1998; Brown et al. 2007). In Section 2, we discuss the exis- 
tence of VMSs in YMSCs and estimate the number of YM- 
SCs formed in the Galactic disk during a certain interval of 
time. In Section 3, we estimate the typical ejection velocity 
produced via exchange encounters between binary stars and 
a VMS. In Section 4, we compare this estimate with the 
results from numerical three-body scattering experiments. 
The discussion is given in Section 5. 



2 VERY MASSIVE STARS 

There is emerging observational evidence that the majority 
(if not all) of massive stars form in a cluster environment 
(e.g. Zinnecker & Yorke 2007 and references therein; cf. de 
Wit et al. 2005; Schilbach & Roser 2008; Gvaramadze & 
Bomans 2008b) and that the mass of the most massive star 
in a cluster is correlated with the mass of the cluster itself 
(Elmegreen 2000; Weidner & Kroupa 2004, 2006). Observa- 
tions also suggest that the maximum mass of ordinary stars 
is saturated at ~ 150 Mq for Md > lO'^ Mq, where Md is 
the mass of the cluster (Weidner & Kroupa 2006), a fact 
which points to the existence of the upper cut-off of stellar 
masses (Weidner & Kroupa 2004; Figer 2005; Oey & Clarke 
2005). 

A somewhat higher upper limit on the maximum stel- 
lar mass follows from the recent work by Yungelson et al. 
(2008), who argue that the birth masses of some of the ob- 
served stars could be as large as ~ 200 Mq (cf. Oey & Clarke 
2005). It is also possible that even more massive stars could 
originate from the coalescence of binaries whose components 
have masses close to the upper limit (Yungelson et al. 2008; 
cf. Leonard 1995). Moreover, numerical simulations of the 
dynamical evolution of YMSCs show that runaway colli- 
sions and mergers of ordinary massive stars can result in 
the formation of VMSs with masses as large as > 1000 Mq 
(e.g. Portegies Zwart et al. 1999; Portegies Zwart & McMil- 
lan 2002; Giirkan, Freitag & Rasio 2004; Freitag, Giirkan & 
Rasio 2006; see also Suzuki et al. 2007). 

Two necessary conditions should be fulfilled for the run- 
away growth of a VMS. First, the parent cluster should be 
compact (e.g. Gaburov, Gualandris & Portegies Zwart 2008; 
Ardi, Baumgardt & Mineshige 2008), i.e. dense enough to 
ensure that stellar collisions are frequent. Second, the clus- 
ter should be massive (~ 10^ — 10^ Mq), i.e. should contain 
a large number of OB stars - the main building blocks for 
the VMSs. Let us discuss these conditions in more detail. 

It is conceivable that the first condition is fulfilled for 
the majority of YMSCs since their characteristic radius at 
birth (i.e. during the embedded phase) is < 1 pc, indepen- 
dently of their mass (e.g. Kroupa & Boily 2002 and refer- 
ences therein). The high central densities in young clusters 
could be either primordial (e.g. Murray & Lin 1996; Clarke 
& Bonnell 2008) or caused by dynamical mass segregation 
(e.g. Portegies Zwart et al. 1999; Giirkan et al. 2004). The 
existing observational data on YMSCs do not allow us to dis- 
criminate between these two possibilities. The observed top- 
heavy initial mass functions in the central parts of two of the 
most massive and dense YMSCs in our Galaxy, NGC 3603 
and Arches, could either reflect the initial conditions in the 
cores of these clusters (e.g. Harayama, Eisenhauer & Mar- 
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tins 2008; Dib, Kim & Shadmehri 2007; see, however, As- 
censo, Alves & Lago 2008) or have a dynamical origin (e.g. 
Harayama et al. 2008; Portegies Zwart et al. 2007). It should 
be noted that the central stellar densities observed in these 
two clusters arc lower than those required for runaway colli- 
sions to occur. It is possible, however, that either NGC 3603 
and Arches were coUisional at birth and their central den- 
sity was already reduced due to dynamical decay or they are 
currently on the way to core collapse and runaway collisions. 
In the following, wc assume that, at a certain stage in the 
early evolution of YMSCs, the star density in their cores was 
high enough for the runaway collisional process to proceed 
and for the formation of VMSs. 

Depending on whether the mass segregation in young 
clusters is primordial or of dynamical origin, the runaway 
collisional process (accompanied by formation of VMSs and 
ejection of high-velocity runaway stars) starts respectively 
at the very beginning of the cluster evolution or after several 
Myr (i.e. when the star density in the cluster's core becomes 
enhanced via the Spitzer instability). In the first case, the 
kinematic age of ejected stars will be close to the age of the 
parent cluster, while in the second one, it will be smaller by 
several Myr. High-precise proper motion measurements for 
OB stars ejected from YMSCs are required to discriminate 
between these two possibilities (e.g. Gvaramadze & Romans 
2008b). 

The fact that our Galaxy is populated by YMSCs (the 
second condition for the formation of VMSs) was acknowl- 
edged only recently (Paper I and references therein; see also 
Figer 2008). The recent discovery of several clusters with 

> 10" M© (Borissova et al. 2006; Figer et al. 2006; 
Davics et al. 2007) and the upward revision of masses of 
the already known stellar systems (Knodlscdcr 2000; Clark 
et al. 2005; Homeier k, Alves 2005; Santos, Bonatto & Bica 
2005; Ascenso et al. 2007a,b; Wolff et al. 2007; Harayama 
et al. 2008) indicate that YMSCs arc much more numerous 
than was hitherto known and that many of them cannot be 
detected due to the heavy obscuration caused by the fore- 
ground and circumcluster dusty material. 

Below, we estimate the number of star clusters with 
Mc\ ^ lO** M0 formed in the Galactic disk during a certain 
interval of time, Ai. Since we are interested in the formation 
of high- velocity early B-type (^ 8 Mq) stars (i.e. the progen- 
itors of the majority of NSs) and high-velocity late B-type 
(^3 — 4 M0) stars (representative of the HVSs), we consider 
two time intervals. At ~ 30 Myr (the lifetime of a single low- 
est mass star producing a NS) and At ~ 350 Myr (the life- 
time of a 3 M0 star). Assuming that the star formation rate 
of ~ 7 - 10 X 10"" Mq yr"^kpc"^, derived by Lada & Lada 
(2003) for embedded star clusters in the solar neighbour- 
hood, is representative for the Galactic disk as a whole (see 
Paper I for more details) , one has that ~ 6.6 - 9.4 x 10® Mq 
and ~ 7.7 - 11.0 x lO'' Mq were formed within a circle of 
10 kpc during, respectively, the first and the second time in- 
terval. Then taking a power-law cluster initial mass function 
with a slope of 2 (e.g. Lada & Lada 2003) and assuming that 
the cluster masses range from 50 to 10® Mq, one has that 
~ 70 - 100 and ~ 800 - 1100 YMSCs with > lO" Mq 
were potentially available to produce the observed popu- 
lation of high-velocity early and late B-type stars, respec- 
tively. In the following we assume that YMSCs of mass of 
~ 10" Mq contain ordinary stars as massive as ~ 100 Mq 



(Weidner & Kroupa 2006), while the more massive clusters 
could harbour the VMSs formed through the merging of or- 
dinary massive stars. 

Similar arguments show that ~ 10^ and lO" star clus- 
ters of mass of ~ 10''* Mq were formed in the Galactic disk 
(within a circle of radius of 10 kpc) during the last 30 and 
350 Myr, respectively. These clusters produce stars as mas- 
sive as 50 Mq (Elmegreen 2000; Weidner & Kroupa 2006) 
and could also contribute to the production of high-velocity 
(^ 200 — 400 kms~^) runaway stars (see Section 3 and Sec- 
tion 4). 



3 HIGH- VELOCITY STARS FROM 

EXCHANGE ENCOUNTERS BETWEEN 
BINARY STARS AND A VERY MASSIVE 
STAR 

Close dynamical encounters between (hard) binary stars and 
a VMS can result in the ejection of one of the stars with a 
high velocity. In this process, one of the binary components 
is replaced by the VMS, while the second one is ejected 
(the so-called exchange encounter). For equal mass binary 
components and zero impact parameter, the typical velocity 
attained by the escaper is ~ 1.8yorb, where V^rb is the orbital 
velocity of the ejected star in the original binary (Hills & 
FuUerton 1980). The ejected star gains its kinetic energy 
at the expense of the increased binding energy of the post- 
encounter (newly formed) binary. 

The production of high-velocity stars in exchange en- 
counters between binary stars and a VMS is similar to the 
process of star ejection in the course of close dynamical en- 
counters between binaries and a supermassive black hole 
(Hills 1988). The main difference is that in the first pro- 
cess the distance of closest approach of the binary to the 
central massive body is limited by the radius of the VMS, 
Ryms, while in the second one it cannot be smaller than 
the tidal radii of the binary components in the gravitational 
field of the black hole. Our goal is to check how this dis- 
tinction affects the typical velocity of ejected stars as well 
as the fraction of encounters resulting in high-velocity ejec- 
tions. Below, we use the results of Hills (1988) to estimate 
the ejection velocity produced in exchange encounters and 
then in Section 4 we compare this estimate with results from 
numerical three-body simulations. 

In the process of a close encounter between a binary 
and a VMS, the VMS can be treated as a point mass if its 
radius is smaller than the binary tidal radius, 

n ~ — a, (1) 

where Mvms, rni and are, respectively, the masses of the 
VMS and the binary components (mi ^ m2), and a is the 
binary semi-major axis. Evolutionary models of VMSs show 
that these stars develop a core-halo configuration (more 
pronounced at the upper end of masses), in which most 
of the mass is concentrated in a dense and compact core 
while the rest of the mass is spread in an extended tenu- 
ous halo (Ishii, Ueno, & Kato 1999; Yungolson ot al. 2008). 
According to Yungelson et al. (2008; also Yungelson, per- 
sonal communication), 99 per cent of the mass of 500 and 
1000 Mq VMSs is confined within a sphere of radius ~ 30 
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and 40 R0, respectively. (In Section 4 we use these figures as 
input parameters for our numerical simulations.) For illus- 
trative purposes, we consider exchange encounters between 
a VMS and a hard binary consisting of main-sequence stars 
with masses mi = 40 M© and m2 = 8 Mq and semi-major 
axis a ~ 3ri ~ 30 R0 (typical of tidal binaries; e.g. Lee & 
Ostriker 1986), where 



ri 



0.8 



mi 



(2) 



is the radius of the primary star (Habets & Heintze 1981). 
For these parameters, one has from equations ((T} and ^ 
that rji"" > 2_RvMS, so that the VMS can be treated as a 
point mass. 

The typical ejection velocity at infinity attained by the 
escapers in the course of exchange encounters is (Ifills 



640 km s ^ a 



100 M, 



1/6 



30 Rp 



nil + m2 



-1/2 



1 1/3 



(40 + 8) Mq 



(3) 



where a is a non-monotonic function of the dimensionless 
closest approach parameter 



I>mi„ = 79 



Rn 



(4) 



and -Rmin is the closest approach distance between the bi- 
nary and the VMS; a ~ 1 for -Dmin — 0, then it reaches 
a maximum of ~ 1.2 for -Dmin — 30 — 40 and then mono- 
tonically decreases with increasing Dmin- For 7?min ~ Rvms 
and a VMS of mass A/vms ~ 500 and 1000 Mq, one has 
from equation (|3]) and ([4]) that Dmin — 40 and Voo is close 
to its maximum value of ~ 1000 and 1100 kms~^, respec- 
tively. Note that the weak dependence of Voo on Mvms im- 
plies that exchange encounters with ordinary stars of mass of 
~ 100 Mq would in principle be sufficient to produce HVSs 
(cf. Section 4). 



4 NUMERICAL EXPERIMENTS 

In this section, we perform numerical simulations of three- 
body encounters in order to obtain the velocity distribu- 
tions for runaway stars produced in the course of exchange 
interactions between hard binary stars and a massive com- 
pact body, either a VMS or an ordinary star of mass of 
50 — 100 Mq. The simulations are carried out using the 
sigmaS package, which is part of the STARLAlfl software 
environment (McMillan & Hut 1996; Portegies Zwart et al. 
2001). For a detailed description of the setup of the scat- 
tering experiments see Gualandris et al. (2005). During the 
simulations we allow for physical collisions when the distance 
between any two stars is smaller than the sum of their radii. 
For ordinary stars we use the mass- radius relationship given 
by equation ([2} , while for the VMSs of mass of 500 Mq and 
1000 Mq we adopt the radii -Rvms ~ 30 Rq and 40 Rq, re- 
spectively (see Section 3) . The initial eccentricity of the stel- 
lar binary is randomly drawn from a thermal distribution in 
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Figure 1. Branching ratio for the outcome of encounters between 
a (40, 8) Mq binary and a single VMS as a function of the VMS 
mass. The different outcomes are; merger (circles), fly- by (trian- 
gles), ionization (stars) and exchange (squares). The error bars 
represent the formal (Icr) Poissonian uncertainty of the measure- 
ment. 



the allowed range (0 — emax), where emax — 1 — 2(ri+r2)/a is 
chosen as to avoid a collision at the first pericentre passage 
(r2 is the radius of the secondary star). The relative velocity 
at infinity between the centre of mass of the binary and the 
central massive body is set to 5 kms~^, in accordance with 
typical dispersion velocities in YMSCs. 



4.1 High- velocity early B-type stars 

First, we focus on exchange encounters producing high- 
velocity early B-type stars. From numerical three-body scat- 
tering experiments, it is known that the least massive binary 
component is more likely to be ejected with high velocity if 
its companion star is much more massive and/or the binary 
semi-major axis is small (e.g. Gualandris & Portegies Zwart 
2007). For illustrative purposes we consider encounters be- 
tween a binary consisting of two main-sequence stars with 
masses mi = 40 Mq and m2 = 8 Mq and a (very) massive 
star with AfvMS ranging from 50 to 1000 Mq. In order to 
maximize the ejection speed, we assume that the binary sys- 
tem is very tight, e.g. formed by tidal capture. In this case, 
a ~ 3ri ~ 0.15 AU (see Section 3 and cf. Paper I). 

In Fig.[T] we present the probability of different out- 
comes (branching ratios) as a function of A/vms- For each 
value of Mvms we perform a total of 10000 scattering ex- 
periments, which result either in a fly-by, an exchange or 
a merger. Ionizations never take place as the binary is too 
hard to be dissociated by the VMS. Mergers occur in a large 
fraction (~65 per cent) of encounters due to the small semi- 
major axis of the binary and the finite size of the stars. Ex- 
change interactions occur in about 6 per cent of encounters. 
During these interactions, one of the binary components is 
captured by the VMS while the other is ejected, sometimes 



On the origin of high-velocity runaway stars 5 



1 — 1 — I I I I I I 




Figure 2. Velocity distributions at infinity for escaping (8 M©) 
stars in encounters between a binary consisting of a primary 
star with mass mi = 40 Mq and a secondary star with 
mass ni2 = 8 Mq, and a single VMS of mass MyMS = 
50 M0,1OO M0,5OO M0 and 1000 Mq (from left to right). The 
binary semi-major axis is a = 0.15 AU. 

with a high velocity. These encounters are the relevant ones 
for the production of high-velocity runaway stars. 

Fig.[2]shows the velocity distribution for 8 Mq escapers. 
As expected, the more massive VMSs are more likely to eject 
stars with high velocities. Fig. [2] also shows that the typical 
ejection velocities attained by the escapers are consistent 
with the predictions derived in Sections. 

In Fig. [3] we show the probability of exchange encoun- 
ters resulting in ejection of the 8 Mq binary component 
with a velocity from 200 to 700 kms~^ (top to bottom). For 
MyMS 100 Mq about 3 per cent of all encounters produce 
runaways with peculiar velocities ^ 400 kms^^. It can be 
seen (see also Fig.[5| that even an ordinary star of mass of 
50 Mq can occasionally (in ~ 1 per cent of encounters) pro- 
duce an escape velocity ^ 400 kms"'^. In order to produce 
escapers with velocities typical of HVSs {Vao ^ 700 kms~^), 
a VMS of several hundred solar masses is required. In the 
case of a 200 - 300 Mq VMS, > 2 per cent of all encounters 
result in an escape velocity of ^ 700 kms~^. This fraction 
gradually increases to > 3 per cent for the more massive 
VMSs. 

Fig. U shows the average recoil velocity, (Vca), of esca- 
pers (8 Mq stars; solid symbols) as a function of the binary 
semi-major axis for four different values of the VMS mass 
MvMS = 50 Mq, 100 Mq, 500 Mq and 1000 Mq. The empty 
symbols indicate the velocity Knax for which 1 per cent of the 
escapers have Voo > V^max- The average and the maximum 
velocities increase with the mass of the VMS, as expected 
from energetic arguments. On the other hand, the fraction 
of high-velocity escapers decreases rapidly with increasing 
a. For MvMS = 50 — 100 Mq, escapers are ejected with ve- 
locities ^ 200 — 400 kms"'^ if the binary semi-major axis is 
in the range 0.15 AU < a < 0.6 AU. The same velocities 



Figure 3. The probability of exchange encounters between a 
(40,8) Mq binary (with a = 0.15 AU) and a single VMS re- 
sulting in ejection of the 8 Mq star with a velocity from 200 to 
700 kms~^ ( top to bottom). 



= 40 Mq 




0.1 1 
a (AU) 

Figure 4. Average recoil velocity of escapers as a function of the 
initial binary semi-major axis in the interaction of a (40, 8) Mq 
binary star with VMSs of different mass: AfyMS = 50 Mq 
(diamonds), AfvMS = 100 Mq (triangles), MyMS = 500 Mq 
(squares), MyMS = 1000 Mq (circles). Solid symbols represent 
the average velocity obtained from a set of 10000 scattering ex- 
periments while the empty symbols indicate the velocity Vmax 
for which 1 per cent of the escapers have Vcc > Knax- The error 
bars indicate the Icr deviation from the mean. For clarity, we only 
show them for one data set. 
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could be achieved with wider binaries (with a up to 1 AU) if 
A^VMS ^ 300 Mq . To produce hypervelocity escapers Mvms 
should be > 500 Mq and the smallest possible semi-major 
axes are needed (a ~ 0.15 — 0.2 AU). Fig. |4] also shows that 
even a 100 M0 star can occasionally produce hypervelocity 
escapers, but the fraction of these events is very small (see 
also Fig.^. 

Note that the average recoil velocity of escapers pro- 
duced in exchange encounters with a VMS is somewhat 
larger than that produced in encounters with an IMBH of 
the same mass (cf. Fig. |4] with Fig. 4 in Paper I). This seem- 
ingly "incorrect" result can be understood if one takes into 
account that the ejection velocity is a non-monotonic func- 
tion of the closest approach of the binary to the central mas- 
sive body (either a VMS or an IMBH) and that for encoun- 
ters with the VMS -Rmin cannot be less than i?vMS, while 
for encounters with the IMBH i?min is limited by the tidal 
radius of the ejected star, ~ (MiMBH/m*)^'^^r*, where 
AfiMBH is the mass of the IMBH and m* and are the 
mass and the radius of the star. For the parameters adopted 
in our simulations, one has from equations ([T]) and ^ that 
Djnin — 40 — 50 and ~ 15 for encounters with the VMSs 
and the IMBHs, respectively. Thus, in the first process the 
smaller Rmin the larger Voo , while in the second one the clos- 
est possible approach to the black hole does not produce the 
maximum ejection velocity (see Section 3; see also Fig. 4 in 
Gualandris & Portegies Zwart 2007)), which in turn results 
in the smaller {Vcx>). 
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Figure 5. Branching ratio for the outcome of encounters between 
a (10, 4) Mq binary and a single VMS as a function of the VMS 
mass. The different outcomes are; merger (circles), fly- by (trian- 
gles), ionization (stars) and exchange (squares). The error bars 
represent the formal (Icr) Poissonian uncertainty of the measure- 
ment. 



4.2 High-velocity late B-type stars 

Now we simulate exchange encounters producing high- 
velocity late B-type stars. In order to derive the probability 
of obtaining the largest possible ejection velocities, we con- 
sider encounters between a very tight binary (mi = 10 Mq, 
7712 ~ 4 M0 and a ~ 3ri ~ 0.06 AU) and a (very) massive 
star of mass of 50 — 1000 Mq . The branching ratio for these 
encounters (Fig. [5} is almost identical to that given in Fig.[T] 
The only difference is the somewhat smaller percentage of 
exchange encounters due to the smaller semi-major axis of 
the binary. 

The velocity distributions for 4 Mq escapers are shown 
in Fig. [S] for four different values of the VMS mass: Mvms = 
50, 100, 500 and 1000 Mq. The typical ejection velocity ob- 
tained for each set of parameters is consistent with the es- 
timates derived from equation ^ . Fig. [7] shows that for all 
values of Mvms about 2 per cent of encounters results in pe- 
culiar velocities ^ 300 — 400 kms~^ and that about the same 
percentage of ejected stars attains velocity ^ 700 kms~^ if 
Mvms > 200 Mq. 

The average recoil velocity of escapers as a function of 
the binary semi-major axis is shown in Fig.|8l As in the case 
of encounters between (40 + 8) Mq binaries and a VMS, the 
average and the maximum velocities increase with the mass 
of the VMS. Both velocities, however, can reach somewhat 
higher values due to the smaller possible semi-major axis of 
the binary (cf with Fig.|3|. Note also that although (Voo) 
decreases with increase of a, the fraction of escapers with a 
given velocity (< (Voo)) could be larger for wider binaries 
(this is because the percentage of encounters resulting in ex- 
changes increases with increase of a; see e.g. Fig. 1 in Gua- 
landris & Portegies Zwart 2007). This efi'ect is illustrated 
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V, (km s-') 

Figure 6. Velocity distributions at infinity for escaping stars in 
encounters between a binary consisting of a primary star with 
mass mi = 10 Mq and a secondary star with mass m2 = 4 Mq, 
and a single VMS of mass Mvms = 50 Mq, Mvms = 100 Mq, 
500 Mq and 1000 Mq (from left to right). The binary semi-major 
axis is a = 0.06 AU. 
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Figure 7. The probability of exchange encounters between a 
(10,4) M0 binary (with a = 0.06 AU) and a single VMS re- 
sulting in ejection of the 4 Mq star with a velocity from 200 to 
700 kms-i (top to bottom). 



(mprrij) + M,,„s - 




a (AU) 

Figure 8. Average recoil velocity of escapers as a function of the 
initial binary semi-major axis in the interaction of a (10, 4) Mq 
binary star with VMSs of different mass: Mvms = 50 Mq 
(diamonds), A/vMS = 100 Mq (triangles), Mvms = 500 Mq 
(squares), MyMS = 1000 Mq (circles). Solid symbols represent 
the average velocity obtained from a set of 10000 scattering ex- 
periments while the empty symbols indicate the velocity Vmax for 
which 1 per cent of the encounters have Voo > V'max. The error 
bars indicate the Icr deviation from the mean. For clarity, we only 
show them for one data set. 



Figure 9. The same as Fig. [7] but for binaries with semi-major 
axis a = 0.1 AU. 

in Fig. [51 which shows that for binaries with a — 0.1 AU 
the percentage of escapers with Voo ^ 400 kms~^ is about 
2.5 times larger as compared to the case of the more tight 
binaries (cf. Fig.!?}. 



5 DISCUSSION 

We performed numerical simulations of dynamical encoun- 
ters between hard massive binaries and a VMS, in or- 
der to explore the hypothesis that this dynamical pro- 
cess could be responsible for the origin of high-velocity 
(> 200 — 400kms~^) runaway stars. In our study we pro- 
ceeded from the similarity between encounters involving a 
VMS and those involving an IMBH (the latter process is al- 
ready known to be able to produce high-velocity runaways; 
Paper I; Gualandris & Portegies Zwart 2007) and the fact 
that the radii of VMSs are smaller than the tidal radii of the 
intruders (the massive binaries), so that the tidal breakup 
and ejection can occur before the binary components merge 
with the VMS (Section 3; see also Gvaramadze 2007). Our 
study was motivated by the recent evolutionary models of 
VMSs (Belkus et al. 2007; Yungelson et al. 2008), which 
suggest that VMSs can lose most of their mass via copious 
winds and leave behind IMBHs with masses of < 70 Mq, 
which are not large enough to contribute significantly to the 
production of high- velocity runaway stars (see Paper I) . We 
therefore explored the possibility that a VMS could produce 
high- velocity escapers (either early or late B-type stars) be- 
fore it finished its life in a supernova and formed a black 
hole. We estimated the typical velocities produced in en- 
counters between very tight massive binaries and VMSs 
(with Mvms > 200 Mq) and found that about 3 - 4 per 
cent of all encounters produce velocities of ^ 400 kms~^, 
while in about 2 per cent of encounters the escapers at- 
tain velocities comparable to those measured for HVSs (i.e. 
> 700 kms"^). We therefore argue that the origin of high- 



8 V.V.Gvaramadze, A.Gualandris and S.Portegies Zwart 



velocity 200 — 400 kms^^) runaway stars and at least 
some HVSs could be associated with dynamical encounters 
between the tightest massive binaries and VMSs formed in 
the cores of YMSCs. In this connection, it is worthy to note 
that the theoretical velocity distribution for HVSs produced 
in the Galactic Centre (i.e. through the dynamical processes 
involving the supermassive black hole) predicts the existence 
of a tail of velocities of up to several thousand km s~ ^ , while 
the peculiar velocities observed for all known (~ 20) HVSs 
do not exceed ^ 800 kms~^ (e.g. Sesana, Haardt & Madau 
2007; Lockmann & Baumgardt 2008). Interestingly, the lat- 
ter figure better agrees with the maximum ejection veloc- 
ity (~ 1000 - 1500 kms~i) produced via dynamical throo- 
and four-body encounters in dense star clusters (sec Sec- 
tion 3 and Paper I) . The future proper motion measurements 
for HVSs with GAIA will reveal what fraction of these ex- 
tremely high-velocity stars originated in the Galactic disk. 

We also simulated dynamical encounters between tight 
massive binaries and single 50 — 100 M0 stars - the most 
massive ordinary stars formed in clusters with — 10^ — 
10" M© (Weidner & Kroupa 2006). We found that from 1 to 
~ 4 per cent of these encounters can produce runaway stars 
with velocities of ^ 300 — 400 kms^^ (typical of the bound 
population of high- velocity halo B-type stars) and occasion- 
ally (in less than 1 per cent of encounters) produce hyper- 
velocity (^ 700 kms~^) late B-type escapers. Note that the 
smaller production rate of high- velocity escapers in this case 
could be compensated by an order of magnitude larger pop- 
ulation of clusters containing 50 Mq stars (see Section 2). 

Our explanation for the origin of high- velocity runaway 
stars requires a very dense stellar environment of the order 
of 10^ — 10^ stars pc~^ (see Paper I). The role of this envi- 
ronment is three-fold. First, it makes possible the runaway 
merging process, resulting in the formation of VMSs. Sec- 
ond, it provides suitable conditions for production of tight 
binaries via the tidal capture process and hardens the exist- 
ing binaries, thereby increasing the probability of energetic 
three-body encounters. Third, it ensures that the three-body 
dynamical encounters are frequent, i.e. the production rate 
of high-velocity escapers could be high. We caution, how- 
ever, that whether or not such high densities exist in the 
cores of star clusters remains unclear to date (see Section 2 
and cf. Paper I). This and numerous uncertainties about the 
initial conditions and early evolution of young stax clusters 
precludes us from making any estimates of the production 
rate of high- velocity escapers (cf. Paper 1). 

In Paper 1, we suggested that some extremely high- 
velocity NSs could be the remnants of hypervelocity massive 
stars ejected via strong dynamical three- or four-body inter- 
actions in the cores of YMSCs (i.e. the origin of these NSs 
should not necessarily be connected with asymmetric su- 
pernova explosions). Our suggestion was based on the fact 
that one of the HVSs known at that time, HE 0437-5439, is 
massive enough (~ 9 Mq; Przybilla et al. 2008a) to explode 
as a type II supernova and to leave behind a high- velocity 
NS. It is important to note that the time of flight of this 
star from the Galactic Centre exceeds the lifetime of the 
star and therefore the extremely high peculiar velocity of 
HE 0437—5439 cannot be attributed to the dynamical pro- 
cesses involving the supermassive black hole in the Galactic 
Centre. The most likely birth place of this HVS is one of 
the YMSCs in the Large Magellanic Cloud (Gualandris & 



Portegies Zwart 2007; cf. Przybilla et al. 2008a; Bonanos 
et al. 2008; see however Perets 2008). Recently, an even 
more massive (~ 11 ± 1 Mq) HVS, HD 271791, was discov- 
ered (Heber et al. 2008; see also Carozzi 1974; Kilkenny & 
Stone 1988) whose birth place cannot be associated with the 
Galactic Centre. HD 271791 is the only HVS with measured 
proper motion and all measurements indicate that this stax 
originated in the periphery of the Galactic disk (Heber et al. 
2008). Thus, the high pecuHar velocity of HD 271791 cannot 
be attributed to the ejection mechanism involving the super- 
massive black hole in the Galactic Centre. Przybilla et al. 
(2008b) proposed that HD 271791 was a member of a mas- 
sive close binary system disrupted in an asymmetric super- 
nova explosion and that the secondary star (HD 271791) was 
released at its orbital velocity (~ 400 kms~^) in the direc- 
tion of Galactic rotation (which allowed them to explain the 
Galactic rest-frame velocity of HD 271791, provided that it is 
on the low end of the observed range 530-920 kms-i). One 
can show, however, that to explain the high space velocity 
of HD 271791 within the framework of the binary-supernova 
scenario, the stellar remnant of the supernova explosion [a 
< 10 Mq black hole, according to Przybilla et al. (2008b)] 
should receive at birth an unrealistically large kick veloc- 
ity of 750 - 1200 kms"^ (Gvaramadze 2009). We therefore 
believe that the more likely origin of the peculiar velocity 
of HD 271791 (and other halo early B-type stars; e.g. Ram- 
speck et al. 2001; Martin 2006) is through the dynamical 
processes discussed in the present paper and in Paper I. 

The same dynamical processes could also be respon- 
sible for the origin of early B-type stars observed in the 
halos of nearby galaxies (e.g. Comeron, Gomez & Torra 
2003). We speculate that these stars can meet and ionize 
the cloudlets of cold gas on their way through the halo 
and suggest that the so-called extraplanar HII regions (e.g. 
Tiillmann et al. 2003) could be the Stromgren zones pro- 
duced by high-velocity runaway OB stars (cf. Gvaramadze 
& Bomans 2008a). Our suggestion could be supported by 
the fact that the Hq,- luminosities of two extraplanar HII re- 
gions in the galaxy NGC 55 (located at ~ 0.8 and 1.5 kpc 
from the galactic plane) are consistent with the possibility 
that the ionizing sources of these objects are single BO or 
09.5 stars (Tiillmann et al. 2003). 

High-velocity early B-type stars ejected at large angles 
to the Galactic plane end their lives in the halo and thereby 
contribute to the population of halo NSs. A possible exam- 
ple of a NS formed in the halo is the high Galactic latitude 
compact X-ray source IRXS J141256.0-I- 792204. This ob- 
ject has a high X-ray to optical flux ratio (> 8700), typical 
of isolated NSs (Rutledge, Fox & Shevchuk 2008). If IRXS 
J141256.0+792204 is indeed an isolated NS, then its distance 
from the Sun is ~ 8.4 kpc, that corresponds to the distance 
« ~ 5.1 kpc from the Galactic plane (Rutledge et al. 2008). 
The surface temperature of IRXS J141256.0-f792204 (in- 
ferred from modelling its spectrum as thermal blackbody) 
suggests that this NS should be < 10^ yr old, provided that 
its cooling follows the standard cooling curves. If one as- 
sumes that the NS was born near the Galactic plane, then 
its peculiar velocity should be > 5000 km s~^, which is too 
high to be realistic. To avoid this problem, Rutledge et al. 
(2008) suggested that either IRXS J141256.0-F 792204 is not 
an isolated NS (i.e. the object is much closer to the Sun) or 
its cooling is non-standard (i.e. the NS is much older). In- 
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stead, we suggest that IRXS J141256. 0+792204 could be 
the remnant of a supernova explosion of an early B-type 
star ejected from a YMSC in the Galactic disk with a veloc- 
ity of z/tt > 170 kms~\ where t* < 30 Myr is the lifetime 
of the star. It is obvious that even in the case of a synimot- 
ric supernova explosion the peculiar velocity of the stellar 
remnant (the NS) will be as large as that of its progenitor 
star (Paper I). We argue therefore that the proper motion of 
IRXS J141256.0+792204 could be as small as ~ 4masyr"^ 
[i.e. ~ 20 times less than that suggested by Rutledge et al. 
(2008)] or even smaller if the NS is near the apex of its tra- 
jectory. In the latter case, the direction of the proper motion 
of the NS could be arbitrary. 
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